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ABSTRACT 

Using the Chandra/High Energy Transmission Grating Spectrometer (hereafter HETGS) we have 
detected Doppler motion of Capella's X-ray emission lines in the 6 - 25A wave-band. The observed 
motion follows the expected orbital motion of Capella's primary. This finding implies that the primary 
G8 III star, not the secondary Gl III star in the Hertzsprung gap, has been the dominant source 
of hot 10 6 8 - 10 7 K plasma at least in the last six years. In addition, the results demonstrate the 
long-term stability of the HETGS and demonstrate small uncertainties of 25 and 33 km s _1 in the 
velocity determination with the HEG and MEG, respectively. 

Subject headings: stars: late-type — binaries: spectroscopic — stars: individual(Capella) — stars: 
evolution — X-rays: stars 
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1. INTRODUCTION 

Capella (a Aurigae = HD 34029, a = 5 /l 16 m 41.4 s , 
S = 45° 59' 53" in Epoch 2000.0) is a well-studied binary 
system consisting mainly of two cool giants: G8 III and 
Gl III (often labeled as Aa and Ab, respectively). The 
two giants have very similar masses of 2.69 and 2.56 M© 
(Hummel et al. 1994) 1 . The cooler, more massive G8 III 
star is a core He-burning (CHeB) clump giant, while the 
hotter secondary Gl III star is a rapidly rotating gi- 
ant crossing the Hertzsprung gap for the first time (Pila- 
chowski & Sowell 1992; Barlow, Fekel, & Scarfe 1993 and 
references therein). Generally both gap and clump stars 
show hot (T ~10 K) coronal emission measure distribu- 
tion, although gap stars are often found more deficient 
in X-rays than clump stars (Ayres et al. 1998). 

In the case of the Capella system, it was not clear how 
much each star contributed to the total X-ray luminosity. 
Linsky et al. (1998) found that the contribution of the 
G8 star to the total Fe XXI A1354 flux (forming at T ~ 
10 7 K) equaled that of the Gl star in the 1995 and 1996 
observations with the Hubble Space Telescope/Goddaid 
High-Resolution Spectrometer, whereas the G8 star con- 
tributed less in other UV lines forming at much lower 
temperatures. Johnson et al. (2002), instead, found in 
1999 with the HST/Sp&ce Telescope Imaging Spectro- 
graph that the G8 star exhibited a negligible contribu- 
tion to the same line, indicating that the 10 7 K plasma 
of the G8 star was variable on a long time scale. In con- 
trast to the finding by Johnson et al., Young et al. (2001) 
found in 2000 with the Far Ultraviolet Spectroscopic Ex- 
plorer that the Fe XVIII A974.85 line (forming at T ~ 
10 6 8 K) originated mainly from the G8 star. Since the 
emission measure distribution of Capella sharply peaked 
at T - 10 6 8 K (Brickhouse et al. 2000; Ness et al. 2003; 
Argiroffi et al. 2003; and references therein), it seemed 
natural to expect that the main contribution to Capella's 
X-ray emission originated from the G8 star. However, it 
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1 It makes little difference in the following work if the mass values 
by Barlow et al.(1993) are quoted, instead. 



was not proven definitively. 

To date, no X-ray telescope can resolve these stars spa- 
tially (separation « 0.06"; c.f., Young et al. 2002 on 
spatially separating the two giants in the UV with the 
-ffST/Faint Object Camera). However, a definitive de- 
termination of the origin of Capella's hard X-rays can be 
achieved by following the Doppler motion of bright emis- 
sion lines in the observed Chandra/HETG spectra in its 
six years' operation. 

2. OBSERVATIONS AND DATA PROCESSING 

Capella has been observed 14 times (see Table 1) in 
the Timed Event (TE) mode with the CTiandra/HETGS 
(Canizares et al. 2005) for purposes of calibration and 
data collection for the Emission Line Project 2 . 

Each dataset was obtained through the Chandra X- 
ray Center data archive and was reprocessed with the 
CI AO tools (version 3.2.1); and ancillary effective area 
and grating line response products were generated with 
the CALDB 3.0.1. The reprocessing included the lat- 
est updates on the ACIS-S chip geometry and the MEG 
grating period 3 . No pixel randomization in the output 
of tg_resolve_events was added and the default CIAO 
extraction settings were applied. 

Since the Capella system is a modestly bright X-ray 
source (photon flux fx ~ 2 cps at th order image in 
the HETG wave-band), its th order image is piled up. 
This leads to poor determination of the centroid peak in 
the th order image and hence a poorer zero-wavelength 
solution. To remedy this, the th order position was re- 
derived based on the intersection of the grating arms and 
the O^-order ACIS frame transfer streak. 

3. LINE ANALYSIS 

To determine accurate Doppler velocities of the bright 
emission lines in the Capella spectrum, we followed this 
procedure: (1) select a region of Capella's spectrum con- 
taining strong emission lines; (2) fit the line position of 

2 http: / /cxc. harvard.edu/elp/ELP.html 

3 http://space.mit.edU/CXC/docs/docs.html#acis_s_geom 
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TABLE 1 

Chandra/HETGS Observations of Capella 







Good Time 




Apparent Velocity 


Corrected R etdifil 


BcLry center 


Sp3iCecrs,ft 




Date Start 


Inter v£il 


PhciSC 


V ■ 

v obs 


3<r confidence 


Velocity 


Correction 


Velocity 


ObsID 


(MJD) 


(ksec ) 






(km/s) 




("Vk km/s") 

V v bary x^ill/ 




fO„„ km/s") 


1099 


51418.329 


14.6 


0.398 


-34.9 


(-19.3 r- 


- -52.7) 


-8.65 


26.0 


0.214 


1235 


51418.512 


14.6 


0.400 


-25.3 


(-7.77 - 


- -41.7) 


0.97 


26.0 


0.272 


1100 


51418.694 


14.6 


0.402 


-10.0 


(+8.94 - 


- -29.5) 


16.4 


26.1 


0.325 


1236 


51418.877 


14.6 


0.404 


-13.2 


(+6.31 r 


- -34.6) 


13.3 


26.1 


0.375 


1101 


51419.060 


14.6 


0.405 


-16.3 


(+4.25 <- 


- -32.2) 


10.2 


26.1 


0.420 


1237 


51419.243 


14.6 


0.407 


-15.0 


(+5.01 <- 


- -32.2) 


11.6 


26.1 


0.451 


1103 


51445.257 


40.5 


0.657 


-6.92 


( + 3.04 r 


- -16.4) 


20.4 


26.9 


0.392 


1318 


51446.560 


26.7 


0.670 


-10.2 


(+4.37 - 


- -24.4) 


16.4 


26.8 


-0.182 


0057 


51606.687 


28.8 


0.209 


67.4 


(+80.1 - 


- +54.7) 


40.4 


-27.3 


0.352 


1010 


51951.516 


29.5 


0.524 


26.2 


(+42.7 r 


- +16.2) 


1.57 


-24.2 


-0.468 


2583 


52393.741 


27.6 


0.775 


54.6 


( + 68.1 r 


- +44.6) 


35.4 


-18.8 


-0.485 


3674 


52909.818 


28.7 


0.736 


-4.34 


( + 9.86 r 


- -17.1) 


21.7 


26.7 


-0.623 


5040 


53258.924 


28.7 


0.092 


7.10 


( + 22.6 r 


- -5.04) 


33.6 


27.3 


-0.785 


5955 


53457.553 


28.7 


0.002 


65.7 


( + 80.3 r 


- +49.1) 


39.2 


-26.5 


0.000 



TABLE 2 

Wavelength ranges selected for Doppler measurement 



No. 


Range A (A) 


Selected Emission lines (\i a b) 






for HEG and MEG 


1 


6.0 - 6.50 


Si XIV AA6.1804, 6.1858 


2 


6.6 - 6.72 


Si XIII AA6.6479, 6.6882 


3 


8.3 - 8.55 


Mg XII AA8.4192, 8.4246 


1 


9.12 - 9.35 


Mg XI AA9. 1687,9.2312,9.3143 






Fe XXI A9.1944 


5 


12.00 - 12.20 


Ne X AA12.132, 12.137 






Fe XVII A12.124 


6 


14.9 - 15.1 


Fe XVII A15.014 






Fe XIX A15.079 






for MEG only 


7 


18.85 - 19.1 


O VIII AA18.967,18.988 


8 


21.4 - 22.2 


O VII AA21. 602, 21. 804,22.098 


9 


24.5 - 25.1 


N VII AA24.779,24.785 



the ±l s * orders of the HEG and MEG spectra separately; 
and then (3) calculate this region's Doppler velocity. The 
average value of these velocities from spectral segments 
and grating orders was taken to be the apparent Doppler 
velocity of Capella for the observation. 

Despite the high spectral resolution obtained with the 
HETGS, almost all of the observed emission lines are 
blended (e.g., all hydrogenic lines are a blend originat- 
ing from two transitions to the ground state, Is 2 Si/2, 
from 2p 2 P3/2 and 2p 2 Pi/2 states). Such an unresolved 
line complex was modeled as a blend in order to make a 
precise measurement of its apparent Doppler shift. For 
the baseline spectral model, a three-temperature APED 
plasma model 4 was chosen to represent the Capella spec- 
trum (similar to the model used in Canizares et al. 2000). 
The three temperatures were fixed (at kT ~10 6 3 , 10 6 - , 
and 10 71 ). The three free parameters, the normaliza- 
tion of line emission and the apparent Doppler velocity 
and width, were then fit by minimizing the C-statistic us- 
ing the Levenberg-Marquardt method as implemented in 

4 APED is a model for collisionally excited thermal plasma 
(Smith et al. 2001). 



ISIS version 1.3.0 (Houck & Denicola 2000). All quoted 
errors in the Doppler velocity term correspond to 3<r con- 
fidence levels (see the fifth and sixth columns in Table 1). 

The regions were selected to include bright emission 
lines; the selected spectral regions and emission lines in- 
cluded in the fitting process are tabulated in Table 2. 
Once all the measurements were done, we took all of the 
velocities per wavelength region per grating, rejected the 
highest and lowest (i.e., possibly anomalous) data points, 
and then took the mean as the apparent velocity of the 
system. The measured Doppler widths were small and 
statistically insignificant (<150 km s _1 ). 

Each measurement of apparent Doppler velocity 
needed to be corrected for the barycentric motion of the 
Earth around the center of mass of the Solar system. The 
barycentric cor rection was made using the algorithm of 
iStumnffl IpSfT) . Although we also applied the correc- 
tion for the motion of the spacecraft, it turned out to be 
insignificant for our analysis. Even at its maximum, the 
scale was less than 1 km s _1 for Capella's location on 
the sky. Barycenter-corrected radial velocities Vbary are 
derived as: 

Vbary — ^o6s ~1~ Chary + C sc 

where V bs is the measured apparent Doppler velocity, 
Chary and C sc are the barycentric correction and the cor- 
rection for the motion of the spacecraft, respectively. 

4. DETECTION OF THE DOPPLER MOTION OF CAPELLA 

Figure 1 shows the apparent radial velocities of Capella 
measured without barycentric correction. The data 
points are seemingly randomly scattered, though close 
comparison with the calculated apparent radial motion 
of the primary Capella Aa (the dashed curve) indicates 
that the observed hard X-ray emission lines of Capella 
trace the primary star's motion. The zero orbital phase 
of the Capella system applied in this work is defined by 
Hummel et al. (1994) as 

P (n) = JL>2447528.45(±0.02)+ 104.022 x n 

where n is an integer. Other orbital parameters are de- 
fined as follows: the orbital inclination i = 137.18°, ec- 
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Fig. 1. — Apparent radial velocities of Capella measured with 
the Chandra/HETGS. The dotted line shows the calculated appar- 
ent radial motion of Capella Aa viewed from Earth (including the 
Barycentric, orbital and systemic motion of Capella Aa). 3cr error 
bars are shown in the plot. 




Fig. 2. — The observed radial motion of Capella vs. orbital 
phase after barycentric correction (see the 7th column in Table 1). 
The measured radial velocity clearly follows the trend of Capella 
Aa (primary). 3<r error bars are shown in the plot. 
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Fig. 3. — Changes in relative integrated line flux of hot (Fe XXI; 
filled circle) and cool (three Fe XVII and Ne X; open circle and 
cross) emission lines. The measured line fluxes are scaled relative 
to the lowest flux observed near phase = 0.4. la error bars are 
shown in the plot. 



centricity e = 0.0, the masses of the primary and com- 
panion stars M p = 2.69 Mq and M s = 2.56 Mq, the 
node of ascension fi(2000.0) = 40.8°, the velocity ampli- 
tudes K\ and K2 for the primary and secondary stars 
are 26.05 and 27.40 km s -1 , and the systemic velocity 5 
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Fig. 4. — Two-velocity-component fit to the Mg XII blend at 
8.42A. At the phase = 0.002 (ObsID 5955), barycenter-corrected 
radial velocities for Capella Aa and Ab are expected to be 55 and 
2 km s —1 , respectively. With the line positions fixed, the ratio of 
the normalization terms for Capella Aa and Ab is then derived to 
be roughly 2:1. 

7 = 29.20 km s" 1 . In Figure 2, the same data points were 
barycenter corrected and mapped into and plotted ver- 
sus orbital phase. The velocity vs. phase diagram makes 
it clearer that the X-ray emission lines detected with the 
Chandra/HETGS closely trace the motion of Capella's 
primary star, as opposed to the secondary. 

5. COMMENTS ON SCIENCE AND CALIBRATION 

In six years of monitoring Capella, the Chan- 
dra/HETGS data reveal that the Doppler shifts of the 
bright X-ray emission lines closely follow the primary star 
of the Capella system. This implies that the G8 III pri- 
mary star has been the dominant X-ray source of 10 7 K 
plasma, at least, in the last six years. 

The identity of the dominant source of hot plasma in 
the Capella system has been somewhat controversial. It 
was long suspected that the G8 giant was responsible for 
the hard X-ray plasma, since the secondary star belongs 
to the class of X-ray deficient stars (late F - early G 
giants; see Ayres et al. 1998). This view is strongly sup- 
ported by a Far Ultraviolet Spectroscopic Explorer ob- 
servation of Capella, which concluded that the far UV 
Fe XVIII A974.85 flux (formed at 10 6 8 K) originated 
solely from the G8 star's corona (Young et al. 2001). 
The evidence pointing to the slower-rotating G8 star 
as the dominating source of coronal emission is partic- 
ularly interesting, considering that the Gl star appears 
to dominate the UV emission formed at lower temper- 
atures (T < 10 6 K in the transition region and chro- 
mosphere; see Woods & Ayres 1995 and Linsky et al. 
1995). On the contrary, the detection of strong Fe XXI 
A1354 (formed at 10 7 K) far-UV emission from the Gl 
star with the HST/STIS by Johnson et al. implies that 
the Gl giant might have dominated the hotter plasma 
(c.f. Linsky et al. 1998). Variability in luminosity is 
another issue for concern, as Capella (notably the G8 
component) is known to vary in the hottest UV emission 
lines (Fe XXI to XXIV; see Brickhouse et al. 2000 and 
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Ayres et al. 2003). In X-rays, however, it is considered to 
be a steady source with no prior detection of flares (e.g., 
see ArgirofH et al. 2003). Is it then reasonable to infer 
that the two giants in Capella are remarkably steady X- 
ray sources? Light curves of Fe and Ne emission lines in 
Figure 3 illustrate that this is not necessarily the case. 
Based on the integrated line fluxes measured in Fe XXI 
A12.28, three Fe XVII AA12. 27+17.05+17.10, and Ne X 
Lya A 12. 13 lines, it appears that the line flux from the 
cooler plasma (i.e., Fe XVII and Ne X lines) changes 
by 20 - 50% over time, while the hottest emission line 
(Fe XXI) is more modulated (by a factor of two) than 
the cooler counterpart. 

Furthermore, the two largest modulations of Fe XXI 
in line flux were observed in the last two observations 
of Capella (MJD = 53258.924 and 53457.553, or phase 
= 0.092 and 0.002, respectively). Examining Figure 2 
closely, we notice that these two barycentric radial ve- 
locities are significantly lower than expected. It is sug- 
gestive that the deviation of the last two data points 
is probably due to an increasing level of contamination 
by the secondary's line emission flux originating from hot 
10 7 K plasma. This scenario is supported by the observed 
characteristics of thermal distribution of plasma in sin- 
gle giants analogous to the components of Capella. The 
clump K0 giant /3 Cet, analogous to Capella Aa, shows 
a thermal distribution sharply peaked at T ~ 10 6 ' 8 K 
(Ayres et al. 1998); on the other hand, the Hertzsprung 
gap GO giant 31 Com, analogous to Capella Ab, is char- 
acterized instead by a temperature distribution steeply 
increasing up to the 10 7 K peak temperature and a large 
amount of plasma at even hotter temperatures (Scelsi et 
al. 2004). Assuming that the Capella components have 
analogous thermal structures, it is then reasonable to as- 
sume that the G8 star is the stronger X-ray source, but 
the Gl star significantly contributes to the hotter emis- 
sion at T > 10 7 K. Figure 4 shows that indeed bright and 
isolated Mg XII Lya doublets can be fit well with two 
line components associated with Capella Aa and Ab. In 
the lower panel of Figure 4, the dashed line centered at 
barycenter-corrected radial Doppler velocity +55 km s _1 
is the Mg XII doublets from Capella Aa (G8 III) and 
the dotted line centered at +2 km s _1 is the same dou- 
blets from Capella Ab (Gl III). The composite of the 
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two models is shown as the solid line in the top panel 
of Figure 4. In this preliminary analysis, the intrinsic 
line widths for both components were set to zero. To 
quantify the velocity perturbation in terms of how much 
flux would be required from the secondary star to offset 
the radial velocity, the normalization terms for each line 
component were allowed to vary. The best fit resulted 
from the line flux ratio (G8/G1) of 1.8 : 1 in the Mg XII 
feature. The weighted mean for the two Doppler veloci- 
ties is (55 km/s x 1.8 + 2 fcm/s)/(1.8 + 1) = 36.1 km/s, 
which is reasonably close to the measured radial velocity 
of 39.2 km s" 1 (see Table 1). This finding indicates that 
the secondary star may have contributed as much as one 
third of the flux in the HETGS wave-band during the 
observations for ObsID 5040 and 5955. This analysis of 
variation in integrated line flux is still very preliminary 
and a further discussion is deferred to a future publica- 
tion. 

Lastly we find that the Chandra/ HET GS appears to 
be a very stable instrument in the long run. The 
improvements in the ACIS-S geometry values and the 
MEG grating period (provided in CALDB 3.0.1, ver- 
sion geomN0005) further enhance its capability to de- 
tect a miniscule dynamic motion of an astrophysical 
source. The differences between the measured and ex- 
pected Doppler velocity of Capella Aa generally lie under 
20 km s , which we shall consider as the systematic un- 
certainty in Doppler velocity determination using the en- 
tire coverage of the HETGS. As for absolute wavelength 
determination (i.e., per emission line), inspection of our 
detailed fits indicates that the r.m.s. of the differences 
derived per emission line are 25 and 33 km s _1 with the 
HEG and MEG, respectively. 
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